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1 . 



Introduction 



For an initially-straight slender column which obeys Hooke's 
law, the buckling load is defined to be the concentric compressive 
load under which the straight configuration becomes neutrally 
stable. The buckling load is a function of geometry, material, 
and end conditions, and it takes the form: 






yn_E _ 4 

A zr 



( 1 . 1 ) 



where P is the buckling load, m is the buckling factor, E is the 
modulus of elasticity of the column material, I is the minimum 
centroidal moment of inertia of the column cross sectional area, 
and A. is the column length. 

Columns with varying cross-sections are often used when the 
prismatic column is inadequate. The varying cross-section 
column is used in many applications in order to decrease weight 
or provide additional strength. Early studies of the effect of 
varying cross-sections on buckling loads were done by Dinnik (1) 
and by Bairstow and Stedman (2). Dinnik presented a table for 
columns which tapered symmetrically about a central prismatic 
section. He based his analytical calculations on the assumption 
of various power laws for the stiffness in the end portions of 
the columns. 

Bairstow and Stedman presented a numerical method of solving 
for buckling loads in columns with varying cross-sections. This 
method was later improved by Ratzersdorf er (3). Wianello (4) 
presented a graphical method of solution for buckling loads, and 
Marley (5) introduced numerical integration into Vianello's 
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process. Newmark (6) introduced a number of refinements to this 
process and summarized the entire procedure. Approximate 
solutions based on the energy method were given by Bleich (7), 
and Boyd (0) conducted both theoretical and experimental studies 
on tapered columns. 

The present study uses a numerical method of solution to 
find elastic buckling loads for columns of varying cross-section. 
Specifically, the columns are tapered, hinge-ended, slender, 
initially-straight, and concentrically loaded. A large family 
of columns is studied, and Dinnik's columns, tapered symmetrically 
about mid-length, are members of this family. Dinnik's data are 
checked and augmented by obtaining additional points and by pro- 
viding four significant figure results. The major portion of 
the results, however, provides buckling loads for columns which 
are not symmetric about midlength. These latter results are 
believed to be new. 

The work can best be described by dividing it into three 
phases: selection of a numerical method of solution; adaptation 

and use of the selected method in solving for the buckling loads; 
and discussion and verification of results. 
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2 . 



Selection of a Numerical Method 



The selection of a numerical method involved the investiga- 
tion of the Bairstow-Stedman , the Wianello, the Rayleigh (9), 
and the Rayleigh-Ritz (IB) methods. Each of these methods is 
designed to find the buckling factor for a varying cross-section 
column (11). 

In the discussion of the four methods, it is assumed that 
I , the moment of inertia at a longitudinal distance x, is known 

A 

at all points. It is further assumed that the boundary conditions 
are those of a hinged-end column: 



y = 0 


at 


x = 0 


(2.1) 


y = D 


at 


x =je 


(2.2) 



Both the Bairstow-Stedman and the Vianello methods utilize 
the differential equation for the deflection curve of a slender 
column under a compressive load and with small deflections: 

£ 1 * ^ +■ Py = 0 (2.3) 

where y is the lateral deflection at the longitudinal distance 
x. It is convenient to modify this equation by a form of 
equation (1.1) : 

o _ ~rr> £ Jz 

r S* 

where I^ is chosen to be the maximum moment of inertia in the 
column. Equation (2.3) then becomes: 



In the Bairstow-Stedman method, a value for m is assumed, 
add after the quantity ml„/l is calculated, the fallowing 

u X 

approximate equations are used in a numerical integration: 



/! ■ yt- * x y't- i 

y> ■ /-• 

y*.* = y< : -i + x yi 



(2.5) 

( 2 . 6 ) 
(2.7) 



dx’ 



The primes indicate differentiation with respect to x, and y 1 ’ 
d 2 

y" = — The column is divided into k increments of equal 
dx^ ' 

length or i =k X . The subscript i indicates the longitudinal 

position and is equal to kx/j? , i.e., when i = □, x = 0 and uhen 

i = k, x =#. The initial entries y and y" are equal to zero in 

o o 

accordance with the first boundary condition and equation (2.6). 

In order to begin the integration an initial value for y^ at x = X/2 
is arbitrarily assumed. Then the equations (2.5), (2.6), and (2.7) 
are successively employed until the integration has proceeded 
along the entire length of the column. If, at the completion of 
integration, the deflection corresponds to the second boundary 
condition, i.e., y^ = 0, the correct value of m was assumed. If 
the end deflection does not meet the boundary condition, new values 



of m must -be assumed u htil' the se cond boundary condition is met. 



In the Vianello method, the initial deflection curve is 

assumed. Then, using equations (2.6), (2.7), and (2.5), respec- 

2 

tively, the curve yk /m is obtained by numerical integration. 

By dividing, point by point, the assumed deflection curve by 

2 

this value and multiplying by the known value k , a value of m 
for each point is obtained. If all values of m agree, the correct 
deflection curve was assumed. If they do not agree, the curve 
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2 

yk /m is assumed to be the deflection curve, and the process is 



repeated until all values of m are the same. Since each assumed 
deflection curve is known to be more accurate than the preceding 
one (6), the V/ianello process will eventually converge to a 
stable value of m. 

The Rayleigh method makes direct use of the fact that a 
column is in a state of neutral equilibrium when it is subjected 
to the buckling load. This can be expressed by the condition 
T = V. T is the work done by external loads as the column is 
displaced from a straight to a bent configuration and \1 is the 
strain energy of bending induced in the column as it is displaced. 
Additionally, T = Pv, where v is the decrease in distance between 
the ends of the column as the column assumes the deflection curve 
y. For small displacements, v is approximated by: 




where M is the bending moment. 



and so 



Since M = Py, 



/ 




Equating T and \l and solving for P, 





7 V </* 



or, 




* 



( 2 . 8 ) 
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The solution for m consists of assuming a deflection curue y(x) 
which satisfies the boundary conditions, equations (2.1) and 
(2.2), and integrating to find l\l and D. It is generally 
impractical to evaluate D analytically; a numerical integration 
procedure is more convenient. An accurate value of m can be 
obtained by applying this procedure to a number of deflection 
curves. Of these, the correct or most accurate deflection curve 
yields the minimum value of m. 

The main advantages of this system are speed and simplicity 
of calculation. The choice of a family of deflection curves 

t 

which can be varied to produce the correct deflection curve is 
difficult, however, and becomes even more so with the intro- 
duction of the non-symmetric column. 

Ritz provided an extension of the Rayleigh method which 
eliminates the previous objection. He employed equation (2.8) 
in devising a method to adjust systematically the parameters in 
an assumed deflection curve to yield a minimum value of m. 

The method, however, employs a comparatively large number of 
separate numerical integrations which would become quite cumber- 
some in solutions for m in non-symmetric columns. 

The investigation of the four methods was not intended to 
be exhaustive. Since each method will-.yield a solution, the 
selection of the method to be used depended only upon the 
relative adaptability of the methods to the present study. The 
l/ianello method with its automatic convergence was judged to be 
most suitable to the present study and is particularly adaptable 
to non-symmetric columns. The procedure outlined by Mewmark 
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and the refinements which he introduced for improved accuracy were 
used in the solutions. These refinements are presented in 
Appendix A. 
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3. 



Solution for the Buckling Factor 

The family of columns studied is an extension of that con- 



sidered by Dinnik. The form is shown in Fig. 1. The moment 
of inertia of c^ss-sectional area at distance x from the left- 
hand end is given by: 



JL - [*/*+ 1// \ n 
T z \ 3/1 +6/4 ) 



O x// * j/f 




$// * x/j ^ ($/Jt+ <*/f) 



4 = C J/f+e/t) (9# ’■«/')* */l * '• o 

f 

Dinnik restricted his study to columns symmetric about midlength, 

l 

i.e., b = e, g = d, c =j?/2. In the present study, non-symmetric 
columns are also included. The shape of an individual column is 
uniquely determined by specifying I^/I^, n > > and • It 

is easily shown that the other parameters are given by: 

j/4 = c// - # / Z / 

J/f = / " j/J - A/J? 

= (- 4jr) j/s 



e/l </// 

where /T = &/&)?* 

It is of interest to note that the postulated variation of 
I is attainable through linear taper for n = 1, 2, 3, and 4. 

Thus: n = 1 corresponds to a column of rectangular cross-section 



2D 





V 



a 



with constant depth and linearly varying widthj n = 2 is approxi 
mated by a built-up column consisting of four angles connected 
by diagonals; n = 3 corresponds to a rectangular cross-section 
oolumn with constant width and linearly varying depth; and n = 4 
is represented by a solid truncated cone or pyramid. 

Since the Uianello method automatically converges, the 
assumed shape of the initial deflection curve is not critical. 
However, the choice deserves some attention because the number 
of iterations required depends upon the accuracy of the starting 
assumption. The boundary conditions for a hinge-ended column 
require zero deflection at the two ends. When I /I is unity, 
the tapered column degenerates into a prismatic column with a 
sinusoidal deflection curve. From preliminary calculations and 
observations of the symmetric columns, it was noted that as 
I^/I^ becomes quite small, the deflection curve becomes approxi- 
mately parabolic. Idith these factors in mind, the initial 
deflection curve is assumed to be: 

y ~ (/' )( / ~ I, / I t ) + l,/l z sinf^yr 1 ') (3.1) 

where u = 1 - 2x /£ . Since the amount of variation in the 
deflection curve caused by the introduction of a non-symmetric 
column was unknown, the variation is ignored in equation (3.1). 

The general solution used to find the critical buckling 
factor can be best explained by an example, Table I. A typical 
non-symmetric column is chosen with n = 2, c/j£ - 0.25, a /£ = 0.0 
and = 0.5. For this example, the length is divided into 

8 segments, i.e., A = 0.125 and k = 8. 
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TABLE I - EXAMPLE OF GENERAL SOLUTION 
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8.14 8.11 7.61 7.67 



The computed values of I /I„ are entered on Line 2. Lines 

X 



3, 4, and 5 are self-explanatory. Line 6 is the assumed value 
of the deflection curve of the column uhen subjected to the 
buckling load. Line 7 was obtained by employing equation (2.6). 



Line 7 is the negative of the quotient of the corresponding 
values on Line 6 by those on Line 2. Line Q is the result of 
altering Line 7 with a refinement by IMewmark. This refinement 
is explained in Appendix A. In a symmetric column, the slope 
( y 1 ) is zero at midlength. For a non-symmetric column, this 
condition is not met, but it is nevertheless used to obtain an 
initial value for Line 9. Using this assumption, the initial 
value on Line 9 at x = 0.45J? is set equal to one half the 
negative of the entry at x = 0.5-^ on Line Q. The entries to 
the left of x = 0.45^ are determined from right to left by 
successive subtraction of the intervening numbers of Line 8. 
Line 10 is determined by employing equation (2.5). 



The initial value is obtained from the boundary condition that 
y = 0 at x = □. The following entries are the cumulative sum 
of the entries of Line 9. The incorrect choice of the position 
of zero slope and the incorrect assumption of a symmetric 
deflection curve lead to a non zero entry at x = i on Line 10. 
The second boundary condition specifies that this value be zero. 




( 2 . 6 ) 




(2.5) 
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Line 11 provides a linear correction to Line 10 enabling the 
boundary condition to be met. The entries of Line 11 are com- 
puted by multiplying the negative of the x = Jl entry of Line 
10 by the respective entries of Line 1. The entries of Line 
12 are the sum of the corresponding values of Lines 10 and 11. 
Finally, the values of m on Line 13 are the quotient of the 
corresponding entries on Line 6 by those on Line 12 multiplied 
by k“~ (64). 

Since it is desired to determine m to four significant 
figures (three decimal places), the criterion for convergence 
is that the maximum and minimum values of m on Line 13 be 
within 0.0005 of each other. If, as in Table I, this has not 
occurred, Line 12 becomes the new assumed deflection curve 
(Line 6) and the process is repeated to convergence. At con- 
vergence, the value of m is taken to be the average of the 
maximum and minimum entries of Line 13. 

This general method of solution was programmed for the 
CDC 1604 high speed digital computer. The factors a// , 

c/ Jl , 1 /l nf and n were varied to cover the usable range 
of columns of this type, and the resulting buckling factors 
are listed in Table II. The factor n varies from 1 to 4 in 
integer steps. For each n, c/j? varies from 0.00 to 0.50 in 
increments of 0.05; it should be noted that values of c/JZ 
from 0.5 to 1.0 are merely repetitions of those from 0.5 to 
0.0. For each c /JZ > I^/I^ varies from 0.05 to 1.00 in incre- 
ments of 0.05 and a /JZ varies from 0.0 to 2c/^ in increments 
of 0.1. It should be noted that a column is not defined when 
a/2 Jl is greater than c/J . 
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870 9.870 



4. Discussion of Results 

In a preliminary investigation the factor k uas varied 
from 20 to 160 in the program solution for a column uhich 
appears to have one of the most stringent requirements for 
correct convergence. It uas assumed that an adequate value 
of k for this column is adequate for the remainder of the 
columns. The column is symmetric, i.e., c/J> = 0.5; and 
n = 4, I^/I^, = ^.10, anc * a /-^ = 0.0. The follouing converged 

values of m uere obtained for each k: 



a) 


k = 


20 


m = 4.81242 


b) 


k = 


40 


m = 4.81206 


c) 


k = 


80 


m = 4.81203 


d) 


k = 


160 


m = 4.81203 



It is evident for the six significant figures shoun, that uhen 
k = 80, the process converges to the proper value. The require- 
ment is for only four significant figures, houever, and k = 40 
safely fulfills this requirement. 

The final results in Table II uere analyzed by: verifying 

that the solutions for prismatic columns, i.e., a// = 1.0 or 

1^/12 = 1.0, agree uith theory; verifying that the results agree 
uith the feu analytical solutions for uhich there is a knoun 
deflection curve; and comparing the results uith those buckling 
factors available from Dinnik's solutions. 

For a uniform slender column uith hinge supported ends, the 

2 

buckling factor m is equal to 7f . The buckling factors in 

2 

Table II agree uith the value of 7f to four significant figures 
for all columns uhere either a// or I-^/I ^ is unity. 
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There are four analytical solutions which are both per- 



tinent to this study and are readily obtainable: 



(1) 


c/f =0.0, 


n=2, 


a /4 =0.0, 




1.0 


(2) 


c/4 =0.5, 


n=2, 


a/i =0.0, 


0SI 1 /I 2‘ 


1.0 


(3) 


c/4 =□.□, 


n=4, 


a// =0.0, 


"“/V 


1.0 


(4) 


c/4 =0.5, 


n=4, 


a/4 =0.0, 


0* Ij/Ij* 


1.0 



These solutions were adapted from similar solutions by 
Bleich (12) and Timoshenko and Gere (13) and are presented in 
Appendix B. A comparison of the buckling factors obtained 
from them with those of identical columns from the numerical 
solution in this study is contained in Table III. The agree- 
ment between the two is excellent. The two non-symmetric 
cases have a maximum deviation between analytical and numerical 
results of B.00QB2. This is quite good since the criterion for 
convergence in the numerical program was that the maximum and 
minimum values be within D.0005 of each other. The maximum 
deviation of D.DBD10 in the two symmetric cases is not as small, 
but is still well within an acceptable limit. 

The buckling factors which are applicable from Dinnik's 
work were obtained from Timoshenko and Gere (13) and are 
listed in Table IV/. Under each of these values is the corres- 
ponding value from this study. There is a small deviation in 
these entries in a large number of cases; n = 1 seems to be the 
prime offender. The maximum difference is B.B6. Since the 
present study agrees so closely with the exact solutions, and 
since the same method of solution is employed at all times, 
Dinnik's results seem to be in error by the amount of deviations 
from the present study. 
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analyt ical 


TABLE III 

CHLCK UE THE CRITICAL 


BUCKLING EACTOR 


N = 2 . u 


C/L= . .0 


A/L=0 . 0 


11/12 
.05 
. 10 
. L 5 
. 20 
.25 


NUMERICAL 
RESUL r S 
2.30237 
3.59827 
4.21168 
4.73361 
5. 19807 


ANALYT ICAL 
RESULTS 
2.80235 
3.59826 
4.21168 
4.73360 
5.19807 


. 36 
.35 
. 4 0 
.45 
. 5 ( 


5.62218 
6.01594 
6.38581 
6. 73624 
7.07045 


5.62218 

6.01594 

6.38581 

6.73624 

7.07045 


.55 
. 6' » 
.65 
.70 
.75 


7.39085 

7.69931 

7.99733 

8.286,9 

8.56657 


7.39085 

7.59932 

7.99733 

8.28608 

8.56656 


. 80 
.85 
. 9C 
.95 


8.83959 

9.10586 

9.36595 

9.62039 


8.83959 

9.10586 

9.36595 

9.62039 



N = 4.0 



C / L = .) 



A / L = 0 . 0 



11/12 
. v' 5 
. 1C 
.15 
. 2 * 
.25 



NUMERICAL 
RESULTS 
2.20692 
3.12135 
3.82248 
A . A 1382 
A. 9348 j 



ANALYTICAL 

RESULTS 

2.20691 

3.12104 

3.82248 

4.41382 

4.9348. 



. 30 
. 3 5 
• A v 
.45 
. 5C 



5.40581 

5.83894 

6.24209 

6.62073 

6.97836 



5.40580 

5.83894 

6.24209 

6.62073 

6.97886 



.56 
.60 
.65 
.70 
. 75 



7.31949 

7.64496 

7.95713 

8.25751 

8.54733 



7.31949 
7.64496 
7.95713 
8.2575 .» 
8.54733 



. 8C 
.85 
.90 
.95 



8. 82764 
9.99933 
9.36313 
9.61971 



8.82764 

9.09933 

9.36313 

9.61970 
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fAtiLE III (CONTINUED) 



N = 2 . 0 


C/L=).5 


A/L=0 • 0 




NUMERICAL 


XnTlytical 


11/12 


RESULTS 


RESULTS 


.05 


4.64476 


4.64466 


.1C 


5.39883 


8 - 


.15 


5.93412 


5.93414 


.20 


6. 36478 


6.36479 


.25 


6.73179 


6.73187 


.30 


7. 05536 


7.05549 


.35 


7. 34687 


7.34688 


.40 


7.61359 


7.61359 


.45 


7.86047 


— 7 . 86045 


.50 


R.091.-1 


8.09098 


.55 


8.30782 


8.30779 


. 60 


8.51291 


8*51288 


.65 


8.70784 


8 . 70781 


.70 


8.99386 


8.89383 


.75 


9.07199 


9.07196 


.80 


9.24308 


9.24305 


. 95 


9.40782 


9.40781 


.90 


9.56683 


9.56682 


.95 


9.72069 


9. 7206 J 



2: 

ii 

• 

o 


C / L = '.5 


A/L=0.0 




NUMERICAL 


ANALYTICAL 


11/12 


RESULTS 


RESULTS 


.05 


3.81652 


3.81648 


.1C 


4.81296 


4. 81219 


.15 


5.49075 


5.49078 


.20 


6.02C32 


6.02034 


.25 


6.46052 


6 . 46061 


.30 


6.84050 


6.84055 


.35 


7.17659 


7.17662 


. 4C 


7.47914 


7.47913 


.45 


7. 755 '4 


7.75503 


.50 


8.00926 


8.00924 


.55 


8.24542 


8 . 2454 j 


.60 


8.46626 


8.46623 


.65 


8.67393 


8.67390 


.7C 


8.87014 


8.87011 


.75 


9.05627 


9.05624 


.80 


9.23345 


9.23342 


.85 


9.40263 


9.40261 


. 9C 


9.56461 


9.56460 


.95 


9. 72016 


9.72007 
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TABLE IV 



COMPARISON WITH DINNIK'S VALUES 







a/e 


h'h 


n 


0 


0.2 


0.4 


0.6 


0.8 




i 


6.48 


7.58 


8.63 


9.46 


9.82 






6.484 


7.600 


8.689 


9.476 


9.819 




2 


5.40 


6.67 


8.08 


9.25 


9.79 






5.399 


6.670 


8.088 


9.251 


9.790 


0.1 


3 


5.01 


6.32 


7.84 


9.14 


9.77 






5.007 


6.303 


7.824 


9.144 


9.776 




4 


4.81 


6.11 


7.68 


9.08 


9.77 






4.812 


6.113 


7.680 


9.083 


9.767 

/ 

j 




1 


7.01 


7.99 


8.91 


9.63 


9.82 






7.009 


7.999 


8.916 


9.555 


9.829 




2 


6.3*7 


7.49 


8.61 


9.44 


9.81 






6.365 


7.487 


8.610 


9.444 


9.814 


0.2 


3 


6.14 


7.31 


8.49 


9.39 


9.81 






6.136 


7.296 


8.488 


9.399 


9.808 




4 


6.02 


7.20 


8.42 


9.38 


9.80 






6.020 


7.197 


8.424 


9.374 


9.805 




1 


7.87 


8.59 


9.19 


9.70 


9.84 






7,870 


8.614 


9.250 


9.668 


9.843 




2 


7.61 


8.42 


9.15 


9.63 


9.84 






7.614 


8.430 


9.150 


9.598 


9.839 


0.4 


3 


7.52 


8.38 


9.12 


9.62 


9.84 






7.524 


8.364 


9.113 


9.621 


9.837 




4 


7.48 


8.33 


9.10 


9.62 


9.84 






7.479 


8.331 


9.095 


9.614 


9.836 


■'"Upper 


values in each 


set of numbers are 


from Dinnik' 


's work; 


these from 


the present work 


are directly belou 


• 





TABLE IV (CONTINUED) 







a/j? 


h'h 


n 


0 


0.2 


0.4 


0.6 


0.8 




i 


Q .61 


9.12 


9.55 


9.76 


9.85 






8.604 


9.100 


9.500 


9.751 


9.854 




2 


8.51 


9.04 


9.48 


9.74 


9.85 






8.513 


9.040 


9.469 


9.740 


9.853 


0.6 


3 


8.50 


9.02 


9.46 


9.74 


9.85 






8.482 


9.020 


9.458 


9.737 


9.852 




4 


8.47 


9.01 


9.45 


9.74 


9.85 






8.466 


9.009 


9.453 


9.735 


9.852 




1 


9.27 


9.54 


9.69 


9.83 


9.86 






9.262 


9.511 


9.701 


9.816 


9.863 




2 


9.24 


9.50 


9.69 


9.82 


9.86 






9.243 


9.499 


9.695 


9.814 


9.862 


0.8 


3 


9.23 


9.50 


9.69 


9.81 


9.86 






9.237 


9.495 


9.693 


9.813 


9.862 




4 


9.23 


9.49 


9.69 


9.81 


9.86 






9.233 


9.493 


9.692 


9.813 


9.862 
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The results, as mentioned previously, cover the usable 
range for columns of the type under study. There might be a 
necessity, however, for obtaining intermediate values within 
this range which are not listed in Table II. Therefore, some 
consideration should be given to the method of interpolation. 
Since the data obtained form curves of second order or higher, 
an error is introduced when employing linear interpolation. 

The amount of error depends upon the position of the needed 
value in Table II. The error increases with increasing n, 
decreasing a /& , and decreasing I^/I^,. 

A comparison of linear interpolation with interpolation 
using Bessel's formula is given in Appendix D. Interpolation 
formulas and functions were obtained from Comrie (14). For a 
column with c /£ =0.5, n = 4, I^/I = 0.23, and a /& = 0.15, it 
was found that there is a Q.11% difference between linear and 
second degree interpolation. 
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APPENDIX A 



NEWMARK REFINEMENTS 

Newmark (6) presents refinements which improve the accuracy 
of the Uianello method. He bases these refinements on the 
assumption that y" is a second degree curve, whereas, in the 
basic method, it is assumed to be linear. He assumes that 
this curve takes the form of a distributed loading curve, Fig. 
2(a), and then solves for the reactions at the end of each 
increment (A). He finds that: 

~bi- (~?% 6 r - s) (a. i) 

~h( 3 % * /or 

Rc*=h( 3s + /Or- <£) 

+ (% * ,0 r * s ) (fl ' 2) 
“ 24 ( 7s * 6r ~ %) 

Equation (A. 2) is used in the present study when the curve 
is smooth. UJhen there is a break in the curve, equation (A.l) 
is used. The breaks in the curve appear where the tapered 
sections meet the prismatic section or each other. The curve 
is assumed to take the shape of Fig. 2(b) and the modification 
is : 

(/4- % +■ 6 r +■ 6r,-j-s 

The reaction R is the concentrated angle change ( A y 1 ). 
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APPENDIX B 



ANALYTICAL SOLUTIONS 



(a) Symmetric Column, n = 2, a// = 0.0, O.B * I /I l.o, c/4 =0.5 
The column is of the form shown in Fig. 3, and: 

A = h( 7 V = 4 f 

The differential equation of the deflection curve is: 

Py = 0 



E Z 4%* + 



1 X chC 

or, substituting for I s 



; c/*<‘ 



■ a </*y z 

1 - y = O 



(B.l) 



where 






*4 



The general solution of equation (B.l) with two arbitrary constants 
A and B is: 

/ = SJ~ [A Sln(h /nf) + Bcos(hJtn%) J <B ' 2) 

where 



h - /<=< z - //*■ 

The boundary conditions for a symmetric hinged end column may be 
taken as: 

5 = &Y> y = ° ^ x = ° 



- 4 ; 
= / 



<7^ ^ = //£ 



— = ^) 

' c/f ° 

Substituting these boundary conditions into equation (B.2) 



results in two homogeneous equations for the constants A and B: 



A 



S /O 



[* ] 






B cos'! A yn[~j~) ~o 









A h +■ B/z = O 



BB 



A 

/ \ 

/ v 




c 

a 

P 

t— » 

o 

O 



o 

4 -» ® 

0) o 

i n 

>>^r 

CO Ctf 
rO 



0 ) 

p 

•H 



ai 



Solutions for A and B exist only if: 

t«n\h J ~ Zh = 0 



which has an infinite number of roots h. The smallest root h-^ 
defines the critical buckling load. This can be written as: 

tajl r = 4 

r Tn (1,/Q 

where 

r - i h, In l l,/lj 



With a given ratio of T 



can be obtained and the root h 



1 



is 



h, 



= 2 y 



=(/ +*h, z ) [ I -(I,/1J±] Z 

(b) Nan-Symmetric Column, n = 2, a/jp =0.0, 0.0-^ 1.0, 

c/j? = 0.0 or 1.0. 

The column is as shown in Fig. 4 and, as for the preceding 
column, the general solution is: 

y - /J~[A (hint,) + B Cos (h In p] 

The boundary conditions are: 

f = fi~) z > y = o «/ x = o 

§ = /. o , y = o x * A 

These conditions yield the two equations: 

A shlh In a//J f ] 1- B cos [h In (J./l.) l ]~0 

8=0 

Therefore , 

S/n [h ~ O 
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T 

C 

o 

z 

0 ) 

u 

p 

bO 

r-l 

&-» 






I 



I 



83 



a//- 0.0, c/( 



or 



A - 



Z Tf 



In a/4) 



and 



^ -[r-d./rtTfe 

(c) Symmetric Column, n = 4, a /£ = c /£ =0.5, 0.0 ^ 1.0. 

The column is of the form known in Fig. 3 and: 



L 



-iM 



The differential equation is: 

- P y - O 



EL 



X c/x' 



or 



cl v * * 

J+ X 



-h Py = O 



substituting = / / X* 

-t Py = 0 



£L\ 


vv 


+ L. 


J+ 1 


jt* 


t 




or 






-P 


z Py 


yt* 




t c/t 



+ § y = o 



(B . 3) 



where 



a PJ 



s = 



£ 4 



Equation (B.3) is a form of Bessel's differential equation and 
has the solution: 



where : 



y = t~ i [/?, T-x (St) + B, Y-. (St)] 

y , rsf i = jsfv(st) 

7 ~i Cdt) S/rctTr- 
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so : 



y = rJ^s 



A, cos (St) + B t s //? (S t) 



and the general solution is: 



y ~ ** [X co Jr) * S t s/n (v) 



(B.4) 



The bogndary conditions are: 



- i, 


y = o 




o 

II 

X 


- J , 


o 

ii 




x = A /z 



Substituting these boundary conditions into Equation (B.4) and 
solving for and B^ yields: 



or 



where 



Then, 



f an (SA\ - ~_£_ 
z « n \Z6J) Jl 



ton @ _ —6 

0 Jl 



/? = £!_= 4_ DEL 

/ ZbJ 2T£ J £ A 



= 4(3*-(I,/J k ) L 

(d) IMon-Symmetric Column, n = 4, a// =0.0, c // = 0.0 or 1.0, 

0.0^ Ij/I^ 1.0. 

The column is of the form shown in Fig. 4 and the general 
solution is the same as for the preceding column: 



/ 



= X 



A z cw(f*) + s/o (^J 



(B.4) 
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The boundary conditions are: 

^ ^ 6 , y = o k - o 

X* = J> y = 0 x - / 

Substituting these boundary conditions into Equation (B.4) and 
solving for and yields: 

s!n (*n) = 0 

which has an infinite number of roots Tnf . The smallest non- 
trivial root 7f defines the critical buckling load. Then, 
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APPENDIX C 
PROGRAM LISTING 
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JOd0499F ♦BRANCH 

PROGRAM NLWMARn 

DIMENSION X(11U)» Y(llu), YY (110) * Y Y Y ( 1 1 O') »RAT (110) »T ( 110) » RAT 10(2 
11 ) ♦ CRLOAD ( 2 1 ♦ 12 ) »ADJ (110) 
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YY ( K/2 ) = - Y Y Y ( K./2+1 ) /2.0 
DO 2u L = 1 » K. 

I F ( L-K./ 2 ) 8,20,10 

YY ( K/2-L) =YY ( K./2-L+1 ) - YYY ( K/2-L+1 ) 



Y ( 1 ) =0 • 0 
DO 21 M = 2»iU 

21 Y ( M ) = Y ( M- 1 ) + Y Y ( i*l- 1 ) 
IF ( Y ( K.+ 1 ) ) 22>2 j»22 

22 DINC=-( Y ( k+1 ) ) /C 
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APPENDIX D 



INTERPOLATION METHODS 

A value af m is desired for c /$ - 0.5, n - 4, I^/I^, = 0.23, 
and a/ A = 0.15. 

(a) Second degree interpolation 



Bessel’s formula is modified to: 




m = m Q + pA^ 


+ Hlaii 

4 


• ( A'.' 


A f / \ 

+ A, ) 


it is assumed that 


i the third differences 


h /l 2 


= 0.15, 


a M = 


0.15 


a// 


m 


/ 

A 
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m = 6.005 + 0. 


5 (.654) 


- .062 ( .099) 




= 0.20, 


a/( = 


0.15 
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608 
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0.2 
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624 


16 
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m = 6.5Q9 + 0. 


5 (.608) 


- .062 (.055) 


h'h 


= 0.25, 


a// = 


0.15 
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. n 
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6.461 


536 




0.1 


6.997 


560 


24 


0.2 


7.557 


562 
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0.3 


8.119 







= 6.406 



= 6.890 



= 6.997 + 0.5 (.560) - .062 (.026) = 7.275 
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h /l z 
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a/i = 


0.15 


a/H 
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A' 
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i — 1 

• 

□ 


7.341 
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7.855 
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-9 


0.3 


8.360 






341 + 0. 


5 (.514) 


- .062 (. 


h n z 


= 0.23, 


a/P = 


0.15 


I l /I 2 
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A" 


0.15 


6.406 


484 




0.20 


6.890 


385 


-99 


0.25 


7.275 


323 


-62 


0.30 


7.598 







m = 6.890 + 0.6 (.305) -.060 (.161) = 7.131 

(b) Linear Interpolation 

I l /I 2 = D * 20 ’ b/J ? = 0,15 
m = 6.589 + .304 = 6.893 

I 1 /I 2 = 0.25, a/ $ = 0.15 
m = 6.997 + 0.6 (.560) = 7.277 

Ij/Ig = 0.23, a/p = 0.15 
m = 6.893 + 0.6 (.384) = 7.123 

‘ If 

(c) Comparison between the two methods 

7.131 - 7.123 n 110/ 
difference = ^ ■ ■■ ■ = 0.11% 
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